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ensemble 










L/a 


A40.24 


1.90 


0.0040 


0.150 


0.190 


24 


A60.24 


1.90 


0.0060 


0.150 


0.190 


24 


A80.24 


1.90 


0.0080 


0.150 


0.190 


24 


A80.24S 


1.90 


0.0080 


0.150 


0.197 


24 


B25.32 


1.95 


0.0025 


0.135 


0.170 


32 


B35.32 


1.95 


0.0035 


0.135 


0.170 


32 


B85.24 


1.95 


0.0085 


0.135 


0.170 


24 



Table 1: The ensembles used in this investigation. The notation of ref. is used for labeling the ensembles. 

1. Introduction 

From experiment it is known that the masses of the nine light pseudo-scalar mesons show an 
interesting pattern. Taking the quark model point of view, the three lightest mesons, the pions, 
contain only the two lightest quark flavours, the up- and down-quaiks. The pion triplet has a mass 
of M„ PS 140 MeV. For the other six, the strange quark contributes also, and hence they are heavier. 
In contrast to what one might expect five of them, the four kaons and the rj meson, have roughly 
equal mass around 500 to 600 MeV, while the last one, the rj' meson, is much heavier, with mass 
of about 1 GeV. On the QCD level, the reason for this pattern is thought to be the breaking of 
the Ua{1) symmetry by quantum effects. The rj' meson is, even in a world with three massless 
quarks, not a Goldstone boson. In this proceeding contribution, we discuss the determination of rj 
and rj' meson masses using twisted mass lattice QCD (tmLQCD) with Nf = 2 + 1 + 1 dynamical 
quark flavours. This will not only allow a study of the dependence of the T] , Tj' masses on the light 
quark mass value, but also an investigation of the charm quark contribution to both of these states. 
Moreover, the rjc meson mass can be studied in principle. For recent lattice studies in A^j = 2 + 1 
flavour QCD see 

2. Lattice Action 

We use gauge configurations as produced by the European Twisted Mass Collaboration (ETMC) 
with Nf = 2+1 + 1 flavours of Wilson twisted mass quarks and Iwasaki gauge action [||, The 
details are described in ref. [^] and the ensembles used in this investigation are summarised in 
table [T[ The twisted mass Dirac operator in the light - i.e. up/down - sector reads ^ 

De=Dw + mo + ipLtJi (2. 1) 

and in the strange/charm sector [^] 

Dh =Dv/ + mQ + ijXa ysT^+jXsT' , (2.2) 

where Dy^ is the Wilson Dirac operator. The value of mo was tuned to its critical value as discussed 
in refs. [^, ^ in order to realise automatic ^(a) improvement at maximal twist |j8|]. Note that the 
bare twisted masses pL^,^ are related to the bare strange and charm quark masses via the relation 

m,,, = AicT ± (Zp/2s)M5 (2.3) 
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with pseudo-scalar and scalar renormalisation constants Zp and Zs- Quark fields in the twisted basis 
are denoted by Xth and in the physical basis by 1//^^/,. They are related via the axial rotations 



XH=e'^^^''l'^if,, Xh = %e'''^^''/\ (2.4) 



With automatic iff (a) improvement being the biggest advantage of tmLQCD at maximal twist, 
the downside is that flavour symmetry is broken at finite values of the lattice spacing. This was 
shown to affect mainly the mass value of the neutral pion mass [g [ig, [11]], however, in the case of 
A^j = 2 + 1 + 1 dynamical quarks, it implies the complication of mixing between strange and charm 
quarks. 



3. Flavour Singlet Pseudo-Scalar Mesons in % = 2 + 1 + 1 tmLQCD 

In order to compute masses of pseudo-scalar flavour singlet mesons we have to include light, 
strange and charm contributions to build the appropriate correlation functions. In the light sector. 



1 



{XuXu XdXd) 



one appropriate operator is given by [ ]12| ] 

■^iWs Vu + VdiYs Vd) ^ V 

In the strange and charm sector, the corresponding operator reads 



(3.1) 



¥c 



m- 



Xc 
,Xs , 



(3.2) 



In practice we need to compute correlation functions of the following interpolating operators 

Zq _ _ 

A.V = {¥si75 ¥s) = {XdYsXc - XJ75Xs)/2 - — {XsXc + XcXs)/2 , 



(3.3) 



i¥ci75 ¥c) = iXsiysXs - Xci75Xc)/2 - — {x^Xc + XcXs)/2 

Zp 



Note that the sum of pseudo-scalar and scalar contributions appears with the ratio of renormal- 
isation factors Z = Zs/Zp, which needs to be taken into account properly. Z has not yet been 
determined for all values of j8 non-perturbatively. 

However, for the mass determination, we can avoid this complication by changing the basis 
and compute the real and positive definite correlation matrix 



'nPkH vPhPi, vPhSi, 

ynsi.e Vsi,Pi, VShShJ 



(3.4) 



with the notation 



Ph = ixdrsXc - X«'75Zv)/2 , Sh = {x,Xc + XcXs)/2 



(3.5) 



and rjxY denoting the corresponding correlation function. Masses can determined by solving the 
generalised eigenvalue problem [1131 IT^ 



<^(0 T]W(?,?o) = AW(?,?o) no) r]W(?,?o). 



(3.6) 
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Taking into account the periodic boundary conditions for a meson, we can determine the effective 
masses by solving 

AW(f + l,fo) f.-'«<"'('+i)+e-'«<"Hr-a+i)) 

for m^"), where n counts the eigenvalues. The state with the lowest mass should correspond to the 
77 and the second state to the 77' meson. 

From the components tJq"' 2 of the eigenvectors, we can reconstruct the physical flavour con- 
tents c^"^^ from 



c 



c. 



(») _ 1 ( 



ZT]r + T7^"0/V2 (3.8) 



with normalisation 



At this point the ratio Z = Z^/Zp is needed again. Assuming for a moment that charm does not 
contribute significantly to the rj and rj' states, one can extract the rj-rj' mixing angle (p from 

cos(0)=cf «cl.'\ sini<l>) = -ci''^ = cP (3.9) 

with ((^)) denoting the rj (rj') state. 

4. Results 



We have computed all contractions needed for building the correlation matrix of eq. (3.4). 



For the connected contributions, we used stochastic time-slice sources (the so called "one-end- 



trick" [15]). For the disconnected contributions, we used stochastic volume sources with complex 
Gaussian noise [|T^]. As discussed in ref. one can estimate the light disconnected contributions 
very efficiently using the identity 

For the heavy sector such a simple relation does not exist, but we can use the so called hopping 
parameter variance reduction, which relies on the same equality as in the mass degenerate two 
flavour case (see ref. [ p3| ] and references therein) 

DJ^^ = B- BHB + B{HBf - B{HBf + Dj^\HBf 

with D/, = {\+HB)A, B=\/A and H the two flavour hopping matrix. We use 24 stochastic volume 
sources per gauge configuration in both the heavy and the light sector. 

We use both local and fuzzed sources to enlarge our correlation matrix by a factor two. In 



addition to the interpolating operator quoted in eqs. (3.1) and (3.2), we also plan to consider the 
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Figure 1: (a) Effective masses in lattice units determined from solving the generalised eigenvalue problem 
with fo/fl = 1 for ensemble B25.32. We show the results extracted from a 3 x 3 matrix, (b) squared flavour 
content of Tj for B25.32. 



/-matrix combination ijoYs, which will increase the correlation matrix by another factor of two. 
The number of gauge configurations investigated per ensemble is in most cases around 1200, and 
for ensemble B25.32 is 1500. Statistical errors are computed using the bootstrap method with 1000 
samples. 

In figure [l] we show the effective masses determined from solving the generahsed eigenvalue 
problem for ensemble B25.32 from a 3 x 3 matrix with local operators only. We kept fo/a = 1 fixed. 
One observes that the ground state is very well determined and it can be extracted from a plateau fit. 
The second state, i.e. the tj', is much more noisy and a mass determination is questionable, at least 
from a 3 X 3 matrix. Enlarging the matrix size significantly reduces the contributions of excited 
states to the lowest states and, due to smaller statistical errors at smaller t values, a determination 
becomes possible. The third state appears to be in the region where one would expect the T]c mass 
value, however, the signal is lost at f/a = 5 already, which makes a reliable determination not 
feasible. 

In figure ^ we show the masses of the tj and r\' mesons for the various ensembles we used 
as a function of the squared pion mass. In addition we show the corresponding physical values. 
The scale was set from /^^ and niji using the results of ref. [||]. It is clear that the r\ meson mass 
can be extracted with high precision, while the r]' meson mass requires a larger correlation matrix, 
which is work in progress. The comparison with the corresponding physical values seems to point 
towards good agreement. 

We also determine the flavour content of the two states as explained above. It turns out that 
the T] has a dominant strange quark content (see right panel of figure [T]), while the tj' is dominated 
by light quarks. For both the charm contribution is rather small, however, for the tj it turns out 
to be significantly non-zero. A preliminary determination of the mixing angle eq. ( |3.9| ) yields a 
very stable value of about 60°. Note that this is the mixing angle to the flavour eigenstates. The 
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Figure 2: Preliminary values for m-q and m^i for two )3 -values in physical units as a function of the squared 
pion mass. Filled symbols represent m-rf, open ones m^/. 

computation of the mixing angle with respect to the T70 and rjg states is in progress. 

4.1 Bare Strange and Charm Quark Mass Dependence 

The results displayed in figure |2| have been obtained using the bare values of jx^y and jXg as used 
for the production of the ensembles. Those values, however, did not lead to the correct values of, 
for e.g., the kaon and D-meson masses, see e.g. ref. Moreover, the physical strange and charm 
quark mass values differ between the A and B ensembles. Hence, figure ^ is not yet conclusive with 
regards to the size of lattice artifacts and extrapolation to the physical point. What we can learn is 
that the light quark mass dependence in both states appears to be rather weak. 

We also have two ensembles A80.24 and A80.24s with different bare values for /x^j 5 and a 
retuned value for K but identical parameters otherwise. Ensemble A80.24s is significantly better 
tuned with respect to the physical kaon mass value [Q]. As seen in figure |2[, our results seem to 
indicate that this change in the bare parameters also has a significant impact on the T] meson mass 
value, while the r]' mass is unaffected within the (rather large) errors. 

5. Summary and Outlook 

We presented a computation of T] and r]' meson masses from Nf = 2 + \ + \ Wilson twisted 
mass lattice QCD. The results we obtained so far are rather encouraging, the t] meson mass can 
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be determined with high precision. Also for the r]' meson mass, we hope to be able to give more 
precise results by increasing the correlation matrix under investigation. 

As it is not so easy to tune bare strange and charm quark masses exactly, we shall in the 
future use a mixed action approach for strange and charm quarks. This will not only avoid the 
complication induced by flavour symmetry breaking in the twisted mass formulation, but it will 
also allow to use the efficient noise reduction techniques in the heavy sector. However, this method 
requires a careful investigation of unitarity breaking effects. 
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